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Abstract. We renormalize the Curci-Ferrari model at two loops in the MS scheme in an 
arbitrary covariant gauge. 
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The Curci-Ferrari model is an extension of Yang-Mills theory where the gluon field, A“, is 
massive, [|], ||. However, as is well known in such a massive non-abelian gauge theory, where 
there is no spontaneous symmetry breaking, the theory is not fully consistent. For instance, it 
does not obey unitarity, [||. 0, |j, though it has been shown to be multiplicatively renormalizable 
to all orders in perturbation theory, 00- !§• Further, as the mass for the gluon is included in 
the Lagrangian in the simplest possible manner, there is the possibility that gauge invariance is 
also broken. However, as was shown in the original work of Curci and Ferrari, 0 |H , if one also 
includes a particular form of mass term for the ghost fields when the gauge is fixed covariantly 
in a non-linear way, the mass term in the Lagrangian in fact preserves BRST symmetry, though 
without restoring unitarity which remains broken. Despite the lack of unitarity the Curci- 
Ferrari model has been of interest recently for a variety of reasons. First, there have been 
several investigations into how a mass gap for the gluon can arise, 0ii- For example, in 
]|, H a gluon and ghost mass term akin to that which occurs in the Curci-Ferrari model has been 
argued to be important in generating a mass gap in ordinary Yang-Mills theories. For instance, 
a non-zero vacuum expectation value for (i( A “) 2 + ac a c a ) emerges under certain assumptions 
where a is the covariant gauge parameter and c a and c a are the respective ghost and anti-ghost 
fields. Such an operator represents the simplest local dimension two operator which can be 
constructed in non-abelian gauge field theories. Also, a similar approach has been considered in 
[?J] where the effective potential of the Landau gauge condensate is constructed. There it was 
shown that at two loops the non-perturbative vacuum favoured a non-zero value for the vacuum 
expectation value. Other studies include |l(J. Second, the model itself has provided a useful 
laboratory for studying gauge theories with a non-linear gauge fixing, 10- Third, massive 
Yang-Mills models have been used to model the phenomenology of the strong interactions. For 


example, such a model was used in [12] to try and understand diffractive scattering. Finally, the 
main motivation for our interest in the Curci-Ferrari model stems from the role of the gluon mass 


as a natural infrared regulator for Feynman integrals which arise in loop calculations, [13 


In the renormalization of QCD at high loop orders the standard method of extracting the 
ultraviolet divergence structure of a Feynman integral with respect to the regularization, such as 
dimensional regularization, is to perform the calculation with massless propagators. However, 
for two and higher loop integrals one has to be careful that spurious infrared infinities are not 
generated. These can occur when one reduces the class of Feynman integrals to vacuum integrals 
by expanding in powers of the external momenta after one or more external momenta have been 
set to zero first. To circumvent this potential problem the technique of infrared rearrangement, 
[15, [lfj, is applied where a temporary mass is added to the appropriate propagators in an in¬ 
frared divergent integral. Whilst this has been a hugely successful technique it suffers from the 
drawback that it is performed by hand and hence limited when considering the large number 
of diagrams which will arise when the loop order increases. Moreover, it was not clear until 
recently how one could develop the method for application in an automatic multiloop computer 
algebra programme. One recent approach to address this problem is that of [jl^, [l8|. There each 
propagator is systematically given an infrared cutoff mass so that the resulting vacuum integrals 
are automatically infrared finite. Provided all the vacuum diagrams are computable at that loop 
order then one can automatize the process. Indeed [17, [l8|] provides a full two loop calculation 
of various anomalous dimensions. Given that the Curci-Ferrari model naturally incorporates a 
mass akin to that introduced for infrared rearrangement and remains renormalizable it seems 
appropriate to consider and develop that model as a useful and alternative tool to compute the 
ultraviolet structure of Yang-Mills theories. This is important since the efficient programmes 


such as Mincer, [19], compute the ultraviolet structure of only massless two-point functions at 
three loops, but for n-point functions with n > 2 one cannot naively nullify all bar two of the 
external momenta, similar to infrared rearrangement, and apply Mincer. This is because this 
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procedure will inevitably give rise to spurious infrared infinities which, in dimensional regular¬ 
ization, cannot be distinguished from the desired ultraviolet divergences. Hence, it would seem 
more appropriate to us to apply the Curci-Ferrari model in these circumstances. However, it 
turns out that currently the model has only been renormalized to one loop in [Q], Moreover, an 
earlier calculation, [14], appears to incorrectly determine the renormalization constants which 
must satisfy the appropriate Slavnov-Taylor identities. Therefore, the purpose of this letter is 
to renormalize the Curci-Ferrari model at two loops in the MS scheme using an arbitrary (non¬ 
linear) covariant gauge fixing term and dimensional regularization. This is necessary given that 
the renormalization constants of the fields, coupling, gauge parameter and masses are required 
for the subsequent renormalization of any Green’s function. As a consequence of our calculations 
we will verify the result of ||] that five renormalization constants are necessary to render the 
model finite as opposed to the three suggested in [i|]. 

The Lagrangian of the Curci-Ferrari model is, pi] , 


L = - -G“G a ^ - — (d^A*) 2 - —A“A a v + d u ^d^c a 

2 

— am 2 c a c a - | f abc A“c b d^c c + C ^Lf eab f ecd c?c b (fc d 


(1) 


where 1 < a < Na with Na the dimension of the adjoint representation, D = d^ + zgT a A“, 
G“ v = d^A* - d u A% - gf abc A b ^Al, g is the coupling constant and rn is the gluon mass. The 
ghost mass is related to the gluon mass by the covariant gauge fixing parameter a. The structure 

constants of the Lie group are f abc and c“ d^ c b = (Ad^c b — (d^c a )c b . Following the usual 
procedures the gluon and ghost propagators are respectively 


- 5 ab 


( k 2 + m 2 ) 


(1 - a)Wk v 
( k 2 + m 2 ){k 2 + am 2 ) 


and 


5 ab 

( k 2 + am 2 ) 


( 2 ) 


The unphysical pole at am 2 in the ghost propagator had led to the hope that it would counteract 
the same pole in the gluon propagator to establish a unitary theory. However, it has been shown, 
[J, |2j. |3j, [J, that this is not the case. To renormalize ( lj) the quantities of the bare Lagrangian, 
denoted by the subscript o, are replaced by the renormalized ones, through, 


A a 0 » = y/ZlA a » , C a Q = \J~Z~ c C a , = VZeC" 

fjo — ^g Q ? m Q — Z rn Tfl , Q^o 'Aa ^A ^ • (3) 


We have assumed initially that there are five independent renormalization constants. In ordi¬ 
nary Yang-Mills one has three renormalization constants since clearly there is no need for Z m 
and the Slavnov-Taylor identities ensure that = 1 in our notation to all orders to preserve 
gauge symmetry. Here we will allow for an independent Z a as in Q . To find each of the renor¬ 
malization constants to two loops the gluon and ghost two-point functions are computed first 
which determine Za-, Z a , Z c and Z m . It is worth pointing out that since the gluon two-point 
function gives Za , Z a and Z m the structure of the ghost mass renormalization is determined 
prior to calculating it. This in fact provides a strong check on our results. In addition we have 
computed the ghost-gluon and triple gluon vertex renormalizations to verify that the correct 
scheme independent two loop /3-function emerges in both the gluon and ghost sectors. The 
renormalization group functions are defined by, (see, for example, pCj]), 


7A (a) 

7 m(a) 


<9 In Za , N <91nZ c . da 

^-air ’ 7 ' (a) = ^ ^ 

31nm cHna 

yn- , 7 a(a) = y;- 

ding ding 


( 4 ) 
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where a = g * 2 /(16ir 2 * * ). From (|3|) these imply 


7 A(a) 
la (a) 
7m(®) 


<9 In Za d\nZ A 

H\ a ) - a - + a la( a J- 


(3(a) 


da 
3 In Z a 


1- 


da 
d In Z rn 
dg, 


- 1a(p) 


da 
1 - 


d In Z n 


a- 


da 


-i 


( 5 ) 


Ordinarily in massless Yang-Mills theory one has the simple relation j a (a) = — 1 a(o)- When 
the condition Z a = 1 is not satisfied the more general relation, (]5]) , emerges. 


To compute the ultraviolet structure of the two and three point functions we have followed 
the strategy of [17, 18] of expanding out the Feynman integrals to the terms involving two or 
three external momenta respectively. Given that the Curci-Ferrari model is renormalizable, 
§ ! 10 , it follows that terms involving more external momenta will not contribute to the 
ultraviolet divergences. Hence one is left with massive vacuum two loop Feynman integrals to 
compute where, unlike the calculation of [p~7|. 18], the masses in the propagators are not all 


the same. Such integrals have been widely studied before and we quote the result of |2l|], for 
instance, for the basic two loop topology. In general we have 


Iki (k 2 + m 2 )(l 2 + m 2 )[(k - l) 2 + m 2 ] 


- (m\ + ml + ml) + Ye + \ 7 )) 

+ ^m 2 lnm 2 + ?n|lnm| + m^lnmlj- + 0 ( 1 ) 


for arbitrary masses mj where f k = f d d k/(2n) d , d = 4 — 2e with e the regularizing parameter 
and 7 is the Euler-Mascheroni constant. The integrals involving different powers of a propagator 
are given by differentiating with respect to the appropriate masses. It is important to realise 
the role of the logarithmic terms. For instance, specifying various values for rnf we have 


Iki (k 2 + m 2 )(l 2 + m 2 )[(k — l) 2 + am 2 } 


and 


Iki 


( k 2 + m 2 )(l 2 + am 2 )[(k — l) 2 + am 2 ] 


- (a + 2)m 2 ^^2 + Ye + “ ln ( 47re 7 )) 

2 

+ ^(a + 2) lnm 2 + alna)— + 0(1) (7) 


- (2a + l)m 2 + Y + ~e. ln ( 4?re 7 )) 

2 

+ ^( 2 a + l)lnm 2 + 2 alna)— + 0 ( 1 ) . 


Therefore, the final renormalization constants could potentially contain In a terms. Moreover, 

in writing two and three point functions in the basic form (j|), one uses partial fractions so that 


1 

(k 2 + m 2 )(k 2 + am 2 ) 


(1 — a)m 2 


(k 2 + am 2 ) (k 2 + m 2 ) _ 


(9) 


which introduces, in addition, powers of 1/(1 —a) which also could appear in the renormalization 
constants. Given that In a and 1/(1 — a) n for n > 1 are singular at a = 0 and a = 1 respectively, 
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it might be expected that at two loops there will be a problem in the Landau and Feynman 
gauges respectively. However, since the original gluon propagator, (|[), can be rewritten as 


- 

k^k u / 1 

1 

. (fc 2 + m 2 ) 

m 2 \(k 2 + am 2 ) 

(k 2 + m 2 )) 


( 10 ) 


there ought to be no problems at a = 1 which provides an internal consistency check on the 
computation. Moreover, the potential a = 0 singularity is avoided by the fact that the basic 
one loop integral has the following e-expansion 


am~ 


Ik (k 2 + am 2 ) 


1 + [1 — 7 — ln(am 2 )]e + 0(e 2 ) 


( 11 ) 


Hence, the counterterms which arise from the one loop diagrams ought to cancel out the ln(a)/e 
poles arising in integrals of the type we have discussed. Indeed in this context we avoid the 
usual approach of subtractions by using the method of ^2|. There one computes the Green’s 
functions in terms of the bare parameters and then rescales them at the end of the calculation in 
terms of the coupling constant expansion of the renormalization constants, (|j), where each term 
of the expansion has already been determined. This strategy is appropriate given that we have 
carried out the calculation automatically using a symbolic manipulation programme written in 
Form, |23|, in order to handle the tedious amount of algebra. The Feynman diagrams for such 
an approach were generated using Qgraf, [^|]. 

Consequently, we find the following MS renormalization constants for the Curci-Ferrari 
model, (D, 


Za — 1 + Ca 


+ C 2 


Z n — 


= 



13 \ 1 


a 


+ 


11 a 


16 16 




Z m — 1 + Ca 


+ C\ 


a 
8 ’ 
a 2 
128 


35 

24 


Z„ = 


11 


1 - ^C A - + C 2 a 
fa € 


53a 1435 

192 + 384 ^ 
121 1 
~2T^2 


a + 

171 
6 e 



+ 0(a 3 ) 



+ 0(a 3 ) 


'fa 2 

35 ^ 

1 

fa 2 

a 

95 \ 

V 16 

32 ) 

e 2 

y 32 

32 

96 ) 


+ 0(a 3 ) 


a 2 11 a 

64 + ~64 


+ 0(a 3 ) 


449 \ 1 
192/ e 


a 2 + 0 (a 3 ) 


( 12 ) 


where f acd f bcd = CA <*> ab • As an additional check on our two loop results we note that the 
0 (l/e 2 ) pole terms of each renormalization constant agrees with the value predicted through 
the renormalization group from the one loop pole. Using (|]), these values for the renormalization 
constants lead to the renormalization group functions 

2 

/ y j 4 (a) = C A (3a — 13)— + C\ ^a 2 + 11a — 59^ — + 0(a 3 ) 
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7« («) 

- C A (3a - 26)^ - C\ (a 2 + 17a - 118) ^ + 0(a 3 ) 


7 c(a) 

2 

= C A (a- 3)7 + C\ (3a 2 - 3a - 95) + 0(a 3 ) 

4 V ' 48 


7 m(a) 

2 

= C A (3a-35)^ + C\ (3a 2 + 33a - 449) ^ + 0(a 3 ) 


/3(a) 

= - yC> 2 - ^C\a 3 + 0(a 4 ) . 

(13) 


The result for the /3-function agrees with the original two loop scheme independent result of 
f25| , [fTj . As the gluon mass in (!lj) corresponds to the coupling of a gauge variant, and therefore 
unphysical, operator its mass anomalous dimension is gauge dependent. Further, the result for 
7 m (o) in the Landau gauge agrees with the calculation of @ as an additional check on our 
computation. It is also worth noting that the explicit sum of the gluon wave function and gauge 
fixing parameter anomalous dimensions gives 


7 A(a) + 7 «(a) 


a 


Ca ~7 + C\ (a + 5) 


16 


+ 0(a 3 ) . 


(14) 


As a final check on our results we note that in the Landau gauge 7 Q (a) = — 7 A (a) as in the 
massless Yang-Mills theory and it is trivial to verify that the Landau gauge renormalization 
group functions coincide precisely with the two loop results of [25, 26, ^7|, 22], In this instance 
the gluon propagator, ( 0 ), takes the simple transverse form 




(k 2 + m 2 ) 

Further, if we define the ghost mass as 


(LV 


~1F 


(15) 


2 2 
m c = am 


(16) 


1 _ 1 

then in general, using Z mc = Z m Z\Z a 2 , 


7m c (a) 


'lm(a) 




d In Z a 
da 


+ ay a (a) 


din Z a 
da 


giving to two loops 


7 me (a) 


—Cau — C\ (18a + 95) — + 0(a 3 ) . 
8 96 


(17) 


(18) 


To conclude we have explicitly constructed all the basic renormalization group functions for 
the Curci-Ferrari model, (I). at two loops in MS. These will be important for using this theory 
to renormalize Green’s functions in Yang-Mills theories where one cannot readily apply current 
automatic multiloop programmes to determine the ultraviolet structure of n-point functions with 
n > 3 and where an infrared regularization is also necessary. Further, given that QCD is more 
appropriate for practical phenomenology it would be interesting to extend our computations 
to include massive quarks in addition to the massive gluons and ghosts. Although that theory 
would remain non-unitary, multiplicative renornralizability ought still to be preserved and hence 
the quark extended Curci-Ferrari model would provide a useful tool for renormalizing full QCD. 
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